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Abstract. Transport and scattering phenomena in open quantum-systems with a continuous energy 
spectrum are conveniently solved using the time-dependent Schrodinger equation. In the time- 
dependent picture, the evolution of an initially localized wave-packet reveals the eigenstates and 
eigenvalues of the system under consideration. We discuss applications of the wave-packet method 
in atomic, molecular, and mesoscopic systems and point out specific advantages of the time- 
dependent approach. In connection with the familiar initial value formulation of classical mechanics, 
an intuitive interpretation of transport emerges. For interacting many-particle systems, we discuss 
the efficient calculation of the self-consistent classical transport in the presence of a magnetic field. 
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INTRODUCTION TO TIME-DEPENDENT SCATTERING 

THEORY 

Potential- scattering theory describes systems with a localized scattering region and a 
potential-free region far away from the scatterer. An incoming plane wave arrives at the 
scatterer and the changes of the phase and the direction of the incoming wave reveal 
properties of the scatterer. While this setup is convenient for the description of collision 
experiments with particle beams, in many other situations matter-waves originate from 
a localized region in space where a potential is present. The latter situation is the topic 
of this course. In this case the formulation and solution of scattering theory is best 
achieved in a time-dependent fashion. Instead of asymptotic potential-free regions in 
space, scattering is recast as an initial value problem, which provides a direct connection 
to the initial value problems familiar from classical mechanics. In quantum mechanics, 
the time-evolution of an initial state is governed by the time-dependent Schrodinger 
equation 

[ihd t -H]\Y(t)}=0. (1) 

For propagating an arbitrary initial state, the basic object of interest is the propagator (or 
Feynman kernel) K, which advances the initial state y(r, t') from initial time t' to a later 
time t: 

/oo 
dr'^Mr'/Mr'/) (2) 
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For Hamiltonians which do not explicitly depend on time, many propagators are avail- 
able in analytic form [1]. In contrast, for explicit time-dependent Hamiltonians only few 
results are known, most of which are related to the existence of Ermakov invariants 



[2, 3, 4, 5]. If the initial state \\f is an eigenstate of the Hamiltonian, the propagation will 
preserve the initial shape of the probability density and the state merely acquires a time- 
dependent phase, oscillating with the eigenenergy of the state. If we propagate any other 
state (which in principle can be decomposed into a superposition of eigenstates) the 
initial density distribution will change shape and move in space. A very instructive ex- 
ample is the Moshinsky shutter [6], where the sudden opening of the shutter gives rise to 
transient effects. An important problem concerns the range of validity of the sudden and 
the adiabatic perturbation theories, which should be applicable for fast or slow changes 
of the Hamiltonian, respectively. For the case of the Moshinsky shutter with adjustable 
shutter opening-time the perturbative and the exact results are given and analyzed by 
Scheitler and Kleber [7]. Transient effects are also the topic of a recent comprehensive 
review article [8]. Besides transient effects, many systems display characteristic revival 
phenomena. Revivals bring back part of the time-evolved wave-packet to its initial posi- 
tion and can lead to complete or partial reconstructions of the initial density distribution 
[9]. Revival phenomena have been analyzed in detail for Rydberg states of the hydro- 
gen atom, and also recently for electronic wave packets in mesoscopic systems [10, 1 1]. 
Here, we first give a brief introduction into the time- and energy-dependent propagator 
and Green function of the Schrodinger equation, and proceed to discuss applications of 
the time-dependent picture in the photodetachment of negative ions. The next section 
shows how to solve the time-dependent Schrodinger equation numerically and how to 
obtain the spectrum of molecular systems. In the last section we model transport through 
semiconductor devices and finally discuss the effects on interactions and boundary con- 
ditions on the transport through nanodevices in a magnetic field. 



WAVE-PACKET EVOLUTION AND GREEN FUNCTION 

The energy-dependent counterpart of the time-dependent kernel is the energy-dependent 
Green function 

1 1 nOO 

G(r,r';E) = lim (r|— — — -|r'> = - / dt ^ E+i ^ t / h K(r,t\r',t'), (3) 
?7^o + E-H + ir) m Jo 

where the limit r\ — > + selects the retarded solution. Only few closed forms of energy- 
dependent Green functions are known and the numerical integration over the highly 
oscillatory propagator fails to give converging results. These difficulties are related to 
the representation of the identity in position space as a ^-distribution 

#(r,0|r',0) = 5(r-r'). (4) 

The Fourier transform of the 5-distribution in position representation is a constant 
function in momentum space and signifies that the energy-dependent Green function 
is the solution of the inhomogeneous stationary Schrodinger equation for all energies 

[E-H]G(r,r';E) = 8(r-r'). (5) 

However, in most physical applications the energy and the momenta are limited to a 
certain range of interest and thus it is possible to replace the ^-distribution for example 



by a localized Gaussian function 5(r), which also represents a Gaussian momentum 
distribution. This step can also be done in a more formal way using coherent states, 
or by switching to the Bargmann representation [12]. The imaginary part of the Green 
function encodes the spectrum of the system. This can be seen from the relation 

Eq. (3) gives the local density of states at the location r — > r' 

-±3[G(r,r;£)] = (r\8(E-H)\r) = [ dA^5(£-£ M;A )|0 njA (r)| 2 , (7) 

where we introduced the complete set of eigenstates <j> n ^ with discrete (n) and continu- 
ous (A) eigenvalues. If we replace the ^-distribution by a normalized state 5, we obtain 
the weighted LDOS 
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r' J r5(r')*G(r',r;E)5(r) 



(S\8(E-H)\S), (8) 



which represents the propagation of an initially localized wave-packet in the time do- 
main, since 



Z[(S\8(E-H)\S)]=Z 

and 
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(9) 



\S(t)) = c-' lHt/h \S(0)}. (10) 
The autocorrelation function 

C(t) = (S(0)\S(t)) (11) 

is the cornerstone of the time-dependent approach to scattering theory. For a given 
Hamiltonian H and initial state 15(0)), the autocorrelation function is obtained either 
numerically or analytically [13]. Depending on the physical system under consideration, 
the initial state \S(0)) has a direct physical interpretation [14, 15, 16], or alternatively is 
conveniently chosen as a vehicle to obtain the correlation function [17]. If the initial state 
is normalized, the autocorrelation function starts with the value C(0) = (5(0) |5(0)) = 1. 
In an open system, where eventually all components of the wave-packet leave the initial 
region covered by |5(0)), the autocorrelation function vanishes in the long-time limit 
lim f ^ 00 C(?) = 0. In a closed system this is not the case and the Fourier analysis of 
C(t) reveals the discrete set of eigenenergies corresponding to the set of eigenstates in 
the system which are represented in the wave-packet. For Hamilton operators which 
are maximally quadratic in positions and momenta, the propagator K is given in closed 
form by Kolsrud [18] based on unitary transformations. In this case the propagator is 
equivalent to a linear canonical transformation of positions and momenta. The position 
space representation of linear canonical transformations has been analyzed in detail by 
Moshinsky and Quesne [19]. In all these cases, the quadratic form of the action preserves 
an initially Gaussian density-profile, and only the phase and width do evolve with time. 
The autocorrelation function and the spectrum of quadratic Hamiltonians cover a wide 
range of physical problems, some are discussed in the next sections. 




FIGURE 1. Photodetachment microscopy in an electric field. The left panel shows the two parabolic 
pathways of the same energy linking the ion source S and the detector. In the experiments by Blondel et. al. 
the detector is placed 0.5 m away from the ion source and the resulting electronic wave-function spreads 
out over several mm. The right panel depicts Einstein's law which relates the photon energy (known 
through the laser frequency) and the electron energy (determined by counting the number of interference 
rings). The intersection point at energy Eq determines the electron affinity. 



PHOTOIONIZATION SPECTRA IN EXTERNAL FIELDS, 
CLOSED ORBIT THEORY 

As a first application, we consider the photodetachment of an electron e~ from a negative 
ion X placed in external electric and magnetic. The ion is then radiated with a laser 
beam of photons y with energy hv. The detachment process is described by 

X- + y^X + e~. (12) 

Conservation of energy relates the photon energy hv, the electron energy E(e~), and the 
binding energy of the electron to the atom, the electron affinity Ea(X~) 

E A (X-)+hv = E(e-), (13) 

which is Einstein's law for the work-function applied to negative ions. If the wave-length 
of the emitted electron is large compared to the size of the emitting object (here the 
negative ion), it is possible to replace the emitter by a point source S(r) = C8 (r — r') and 
only consider the orbital characteristic explicitly [20, 15]. We obtain an inhomogeneous 
Schrodinger equation for the electronic wave-function originating at r' 

[E-H] ¥ (r;r',E)=S(r), (14) 

where C denotes the strength of the point source. The photodetachment process is mod- 
elled as a two-step process: first the energy gained by the photon- absorption promotes 
the electron from the bound state to an unbound state. In the second step, the residual 



effect of the remaining neutral atom X is neglected and we consider the propagation in 
the potential given by the external fields, contained in the Hamiltonian H. 

The measurement of the onset of the photodetachment current as function of photon- 
energy has been used to extract the electron affinity, but contains rather large experimen- 
tal uncertainties. A much more precise method takes advantage of interference effects 
in order to measure the energy of the detached electron. Interference effects require the 
presence of external fields in order to construct multiple pathways from the ion to the 
electron detector, along which the electron travels coherently. The simplest configuration 
is the application of a uniform electric field in the Hamiltonian, 

H md =^--eS-r (15) 
2m 

which creates a "virtual double-slit" (Fig. 1). This interferometric setup was proposed 
by Demkov, Kondratovich, and Ostrovskii [21] and experimentally realized by Blondel, 
Delsart, and Dulieu [22]. The observed quantity on the detector is the spatially-resolved 
rate of incoming electrons, the current density j(r,2s) for different laser frequencies and 
thus for different electron energies E. 

The solution of Eq. (14) for a point-emitter is given by the Green function 

y(r;E) = Jdr' G(r,r , ;E)S(r > ), (16) 

and the probability current-density becomes 

j(r; J E) = -3{v/(r; J E)*Vr; J E)}-^^| V /(r; J E)| 2 , (17) 
m m 

— * 

where srf denotes the vector potential. Applying the equation of continuity 

V.j(r;£) = -?3{S(r)>(r;£)}, (18) 

allows us to obtain the total current by integrating over a surface enclosing the source 
Six) 

J(E) = ~Z(S\G\S). (19) 

For a point-like source, the total current is proportional to the local density of states 
given by the imaginary part of the Green function 

J(E) = -hG(r',r';E). (20) 
n 

For the field-free case we have 

v '—2mE 
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G bee (r,r';E) = — — T -exp 
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FIGURE 2. Photodetachment currents for different electric and magnetic field configurations. The 
values of the magnetic field are 38 = 0.5 T and for the electric fields S — 200 V/m, the respective Hamilton 
operators are given in the text. 



and the current near the detachment threshold has a square-root dependence on the laser 
energy 

J(hv) oc®(hv-\E A \) Vhv-\E A \. (22) 

The energy-dependent Green-function is known in closed form for very few problems 
[23], whereas the time-dependent propagator is available for a much larger class of 
problems [1]. For all combinations of uniform electric and magnetic fields, a quadratic 
Hamilton operators results and therefore for initially Gaussian states the autocorrelation 
function can be obtained as analytical expression. Thus one can compute the spectrum at 
least numerically with very high precision for all field configurations. The resulting pho- 
todetachment currents are shown in Fig. 2 and display the strong influence of external 
fields on the photodetachment process, which causes large oscillations around the field- 
free current (a), Eq. (22). The presence of the electric field (eq. 15) causes a stepwise 
modulation of the current (b), which is linked to the appearance of a new interference 
fringe in the spatially resolved current density [14]. The case of parallel electric and 
magnetic fields (e) with Hamiltonian 

p2 1 e gg 

Hg\\3B= 2^-e^z + -mcof(x 2 +y 2 )-COi(yp x -xp y ), (Oi = — (23) 

is discussed in detail in Refs. [24, 25]. For crossed electric and magnetic fields (d) with 

H E ±S8 = -z e^x+-m(of(x 2 +y 2 ) - (Oi{yp x -xp y ), ft)/ = — (24) 



an interesting substructure within magnetic Landau levels emerges, which is analyzed 
in Ref. [26]. 
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FIGURE 3. Photodetachment current of negative sulfur in the presence of a magnetic field. The Zeeman 
splitting of the levels and the thermal occupation number are taken into account and require to superimpose 
the photodetachment currents of several levels (adapted from Ref. [30]) 



An elegant interpretation of the very prominent interference phenomena encoded in 
the total current is given by closed-orbit theory [27, 28, 29, 25]. The basic idea is to 
approximate the Green function, given by the Laplace transform of the time-evolution 
operator 



G(r,r';£) = -- / dt^'^ (r| exp[-iHt/h] \r'} 
n Jo v v ' 

Feynman path integral 

= --J dtc iE '/ h a(ry ,t)cxp[iS(r,t\r' ,0) /h] : 



action 



in the semiclassical limit h — > 0. The saddle points are given by the following condition 
for the action 

^[Et + S(r,t\r',0)] t=tk = 0, (25) 

which selects all classical trajectories corresponding to a fixed energy E. The travel times 
are denoted by t k . The semiclassical Green function reads 

G(r,r ;E) « G(r,r ;E)+— > a(r,r ,t k ) ,,,,,, > 

lh ti V\S(rj k \r',0)\/(2nh) 

(26) 

and represents a sum over classical paths weighted with complex phases. The first term 
represents the contribution of the t = pole of the propagator and requires a careful 
contour evaluation. The total current is related to all orbits returning to the point of 
origin and thus requires to identify the corresponding closed orbits in position space 
(not necessarily closed in momentum space). If besides the first derivative of the action 
also the second derivative vanishes, the primitive semiclassical approximation diverges 
(diffraction catastrophe) and higher order terms are required. The challenging evaluation 
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FIGURE 4. Contour plots of the potential energy surface of CH^I for the ground state (left panel) and 
the excited state (right panel). Countour labels denote the energies in Hartrees. 



of the corresponding diffraction integrals is a topic of catastrophe theory [25, 31]. 
Measurement in parallel and perpendicular fields by Yukich [30] and Blondel et al [32] 
are in excellent agreement with the theoretical description given above. In Fig. 3 we 
compare the theoretically calculated photodetachment current of the sulfur-ion in the 
presence of crossed electric and magnetic fields {£$ « 1 Tesla) with the experimental 
data [30]. 

For coherent quantum sources, which have an extension of the order of the wavelength 
of the emitted particles, the source structure does influence the cross-section and has to 
be taken into account. An example of such a system is the atom-laser from macroscopic 
Bose-Einstein condensates, where a weak perturbation by a radio-field causes a coherent 
outcoupling of atoms from a trapped condensate [33, 34]. Depending on the size of the 
BEC interference phenomena will occur or a tunneling regime prevails with only one 
possible trajectory [14, 15, 16]. 

MOLECULAR PHYSICS 

The time-dependent approach to molecular physics was explored in the seventies by 
Heller [13, 35] and is the topic of a recent textbook by Tannor [36]. Wave-packets 
are extremely valuable tools in theoretical chemistry and have been used to describe 
photochemical reactions, such as the photodissociation of molecules [37]. In quantum- 
chemistry, the potentials are generally not available in analytic form and numerical 
methods are required to propagate any initial state. Moreover, often the molecular kinetic 
energy operator is coupled to spatial coordinates which further complicates the analysis. 

As a specific example we discuss the photodissociation of a linear, triatomic molecule 
following the time-dependent approach of Ref. [37]. We describe the molecular Hamil- 
tonian within the Born-Oppenheimer approximation, which yields the potential energy 
surface (PES) for the motion of the nuclei, given by the addition of the nuclear Coulomb 




time [fs] energy [meV] 

FIGURE 5. Autocorrelation function (left panel) and its Fourier-transform, which reveals eigenstates 
with eigenenergies at Eq = 109 meV, E\ = 176 meV, E2 = 242 meV, £3 = 260 meV, £4 = 326 meV. 



energies and the electronic energy of a specific electronic state for a fixed nuclear con- 
figuration. The electrons are assumed to follow the nuclear motion instantaneously and 
cross-couplings between the electronic and nuclear momenta are neglected. The intro- 
duction of Jacobi coordinates facilitates the discussion and separates the center-of-mass 
motion from the relative motion. To be specific, let as consider methyl-iodine CH3I, 
which has a tetrahedral structure with the 3 hydrogen atoms arranged in a plane, the car- 
bon atom situated above with the iodine attached and pointing outwards normal to the 
hydrogen plane. For simplicity, we will collapse the three hydrogen atoms to a single 
"atom" (denoted by X = 3H). The kinetic energy of the system is then given by the sum 
of the kinetic energies of the three nuclei 

222 
rr _ Pi + _P^ + _P]^_ (27) 



2Mj 2M C 2M X 

Jacobi coordinates for N particles located at positions ri , . . . , r# are specified by 

li = ^-r 2 = n-r 2 , (28) 
t Miri+M 2 r 2 

6 = M, + M 2 <29) 

= M l + ... +Mj (30) 

Together with the total mass and the reduced mass 

N / i \- 1 

m = £m ( -, ^r 1 = E m j) + M i+v ( 32 ) 
i=\ \j=\ / 
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FIGURE 6. Probability densities of the lowest five energy-eigenstates corresponding to the labelled 
peaks in the spectrum in Fig. 5, obtained by propagating an Gaussian wavepacket and projecting out 
eigenstates with Eq. (36). 



the Hamiltonian is expressed as 

h 2 N ~ l h 2 - 
H = Hcn + Hra = "2^V! + £ -— V| + V(ft, . . . (33) 

If we consider CH3I as a linear triatomic molecule, where the postions of the three 
atoms are described by the three coordinates r/,rc,rx, the Jacobi coordinates reduce 
the problem to a two-dimensional one, where only the distances t,\ = yq — rx and 
£2 = r C m(CX) - ri appear. For the atomic weights M(X) =M\= 3u, M(C) =M 2 = 12w, 
M(I) =M 3 = 127m, we obtain & = r cx and £ 2 = r a + j^r cx , Ml = j^u, /l 2 = j^mu. 
For the potential V(4i,^2) we use the expressions given in Ref. [38] 

V gr (r cx ,r c/ ) = -D e -£*+D e [e-^^-^)-l] 2 

+ + (Jtf _jk)e-«(^-^)][r CT -rf. x e-°( rCT -'cr)] 2 , (34) 

with all length measured in Bohr radii as = 5.29 • 10~ n m and the parameters D e = 
0.0874 Hartree, E* = 0.0346 Hartree, a = 0.4914/og, /3 = 0.899/a fi , r e CI = 4.043 a B , 
r cx = 0-6197 as, k = 0.0363 Hartree/a|, k e = 0.1463 Hartree/a|. The potential energy 
surface is shown in Fig. 4. In the first step we propagate an arbitrary wavepacket (for 
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FIGURE 7. Autocorrelation function (left panel) and its Fourier-transform resulting from evolving the 
ground-state y/£ (| ) on the excited state potential energy surface. The continuous spectrum is directly 
proportional to the photodissociation cross-section. 

example a Gaussian of width a = 5 ■ 10 i m) over 4000 timesteps At = 2 ■ 10~ 16 s 

on a grid of [—1.5, +1.5] x [—1.5, +1.5] (10~ 10 m) 2 . Two algorithms are commonly 
used, the split-operator method based on Trotter's formula [39, 40] and the direct 
polynomial expansions of the time-evolution operator or lHt l h [41]. The split-operator 
method relies on the Fast Fourier Transfrom (FFT) method to apply momentum and 
position-dependent operators by simple multiplications to the wavefunction in position 
and momentum representation. The symmetrized split-operator algorithm is accurate up 
to order A? 3 : 

x e -iAf/ftV(r) jp-l e -iAf/»r(p) j? & -iAt/HV{r)] N 

xe iA */ fty W/V(r/)- (35) 

Here, & denotes the Fourier transforms and J^ _1 the inverse Fourier transform. The 
resulting autocorrelation function is displayed in the left panel of Fig. 5. The peaks 
of the Fourier transform of C(t) indicate the postion of eigenfunctions (right panel in 
Fig. 5). At these energies, we project out in a second run the eigenstates by recording 

T 

y Ei {i)= I dfe-^/VCl,')- (36) 
Jo 

Fig. 6 shows the lowest five energy eigenstates. To describe a photoionization process, 
we need the potential energy surface of the excited molecule, where the iodine is in a 
new electronic state 

CH 3 I + r ^CH 3 +I*( 2 P l/2 ), (37) 

separated from the ground state energy by the energy of the photon. The PES of the 
excited state from Ref. [38] reads in Hartrees 

1 r 2 

V a (rcx,r CI ) = 9.6\%&- lA{rci+Q2rcx ^ aB + 2.604c- 1 - 2rc '/ flB + -0.0362-^, (38) 
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FIGURE 8. AlGaAs/GaAs heterostructure with different layers of semiconductors. The quantum well 
forms at the interface between AlGasAs and GaAs and is filled with electrons donated from silicon atoms. 



and is shown in Fig. 4. In addition, also the knowledge of the dipole moments fi^ is 
required, here we choose for simplicity = 1 . The cross-section of the photoreaction 
is given by propagating the ground state y/£ on the excited potential energy surface, 

S(l,t) = e-^+W* [KiYfeb(I)] (39) 
recording the autocorrelation function 

C(t) = |df S*(|;0)S(|;f), (40) 

and finally applying the Fourier transform to C(t), which results in the spectrum dis- 
played in Fig. 7. The excited PES does not support bound states and thus the autocor- 
relation function decays with time lim f ^ooC(?) = 0. This decay implies a continuous 
spectrum for the process under consideration. The molecule CH$I still serves as a pro- 
totype for calculating PES and propagating wave-packets and recent results are given in 
Refs. [42, 43]. 

We have described the complete reaction dynamics using only time-dependent meth- 
ods and without need to diagonalize matrices. The last point is important for the appli- 
cation of the time-dependent theory to transport in mesoscopic systems. 



TRANSPORT THROUGH MESOSCOPIC SYSTEMS 

In this section we consider many-electron systems in semiconductors at low tempera- 
tures. An introduction to the physics of mesoscopic systems present in nanostructures 
is given by Davies and Ferry [44, 45]. AlGaAs/GaAs heterostructures are nanodevices 
consisting of layers of different semiconductor crystals, which have slightly different 
band structures. At the two-dimensional interface between two materials, the mismatch 
in the conduction- and valence-band results in the formation of a potential well, which 
traps and confines electrons along the two-dimensional interface (see Fig. 8). The elec- 
trons in the well originate from implanted donor atoms (for example silicon atoms), 



which are spatially separated from the interface. This separation in combination with 
extremely clean crystal growth using molecular beam epitaxy (MBE) leads to a high 
electron mobility and the suppression of scattering events in the well. At liquid helium 
temperatures (4 K), the electrons can transverse distances of several /im without loss of 
coherence. 

The strong confinement of the mobile electrons at the interface reduces the intially 
three-dimensional problem to an effectively two-dimensional one, which is very well 
suited for using wave-packet techniques. Refs. [17, 46] contain a detailed description 
of the time-dependent approach to transport in mesoscopic systems. Again, the key 
quantity is the time-dependent correlation function. In experiments, magnetic fields are 
commonly used to alter the electronic pathways and to exploit additional magnetic- 
flux dependent phase effects. In the presence of a magnetic field perpendicular to the 
interface, the split-operator method has to be modified since the Hamiltonian contains 
products of the momentum and position operator which require to divide the Fourier- 
transform step into two parts. The vector potential for a homogeneous magnetic field 
SB = (0,0, i^) in the symmetric gauge becomes si = (— y,x, 0)&/2. The Hamiltonian 
splits into three parts 

H = ^ l - a) ^y+2t + co LPyX+-Col(x 2 +y 2 )+V(x,y), G> L =^, (41) 

T px,y t pv,x V eS (x,y) 

where the mixed momentum-position representation for the kinetic energy is possible 
since [p x ,y] = [py,x] = 0. The new propagation algorithm reads 

\l/{r,t' + NAt) « e - iAt / hv ^/ 2 

r i N 

x^^^r,/), (42) 
where & y denote partial Fourier transforms with respect to only one-dimension. 



Magnetic focussing 

The effects of a magnetic field on electron transport are particularly strong in the 
magnetic-focussing configuration shown in Fig. 9. The electric current originates from 
the lower left chamber and has to pass through a Quantum Point Contact (QPC), which is 
in fact not a pointwise object, but rather a narrow constriction. Afterwards the electrons 
are following a circular trajectory, whose radius is proportional to the inverse of the 
magnetic field 

V2mE 

r c = — 7t>~- (43) 

If multiples of the radius r c equal the spacing to another constriction located at x = 
—500 nm, current can enter the left chamber. The macroscopic current flow from the 




FIGURE 9. Left panel: two-dimensional potential landscape of the device. The T-shaped potential 
of two quantum point contacts (at x = —500 nm and x = nm) forms two constrictions, while the 
background charges produce an irregular potential landscape. Right panel: gray-scale representation of 
the probability density and flux (arrows), obtained by starting a wave-packet located at (x,y) = (0,0) and 
Fourier-analyzing the time-evolution at the energy E = 10 meV. The obstacle at (—350, 125) nm interrupts 
the electron flow from the right chamber to the left one and results in a reduction of the transmission. 



right chamber to the left one is therefore strongly magnetic field dependent. So far we 
have neglected any additional potential perturbation caused from residual donor charges. 
In a real device, these charges lead to potential fluctuations. The fluctuations are on a 
scale of 10 percent of the Fermi-energy Ep of the electron system (typical values of the 
Fermi energy are Ep = 10 meV). The effective potential V(x,y) of the device is shown 
in the left panel of Fig. 9. The transport paths at the Fermi-energy through the device are 
revealed by propagating a wave-packet starting in the lower right chamber and tracking 
its time-evolution. By Fourier-transforming the time-evolved wave-function, we obtain 
the stationary state connecting both chambers at the Fermi-energy of the electron gas, 
visible in the right panel of Fig. 9. 

A direct imaging and verification of the theoretically calculated transport-pathways 
is possible using scanning-probe microscopy. In these experiments, the electronic path- 
ways are locally perturbed by inducing a bump in the potential with a metallic tip, placed 
closely above the device. This tip can be moved with an accuracy of 10 nm and the sys- 
tematic mapping of the change of the macroscopic electron flow through the device due 
to the position of the tip allows one to visualize the transport pathways. The experimen- 
tal data of Aidala et al. for a device of area 4/im 2 is shown in Ref. [47], together with the 
theoretical simulations done by wave-packet runs, described in more detail in Ref. [17]. 

Aharonov-Bohm interferometer 

As the second application of the wave-packet method in mesoscopic nanodevices, we 
consider an electron-interferometer in the presence of a magnetic field. The Aharonov- 
Bohm effect [48] allows one to detect the interference of coherent electrons in a two- 
path interferometer as oscillations in the probability density. The oscillation-period 
of the magnetic flux is given by the flux quantum h/e. In the nanodevice version 
discussed here [49], the magnetic field also penetrates the arms of the interferometer, 
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FIGURE 10. Left panel: sketch of the half-circular Aharonov-Bohm interferometer with four attached 
leads. The right panel shows a close-up of the potential at the crossing, where reflections happen. The 
steep potential profile along the arms of the interferometer (lower right panel) supports several populated 
transverse modes below the Fermi energy (£> = 8 — 15 meV). Adapted from Kreisbeck et al. [51]. 



and thus the eigenstates do not only acquire a phase, but are also directly distorted by 
the magnetic field. Using special etching techniques, interferometers with arm-lengths 
in the /im-range can be fabricated [50]. The device in Fig. 10 was designed to facilitate 
the comparison of theoretical and experimental results and special care was taken to 
attach to the half-circle leads with minimal imperfections. If the lengths of the two 
paths si and si differ, an additional wave-number-dependent phase occurs, given by 
Aa = &f(s2 — si), where kf = \J2mEp / h denotes the Fermi wave-number. Ideally, the 
transmission probability along the paths becomes T cos(e(j)/h + Aa), with magnetic 
flux = area ■ SB through the enclosed area. The wave-number can be controlled by a 
perpendicular electric field applied via a top-gate electrode. The simple linear relation 
between wave-number and phase Aa does not take into account time-reversal symmetry, 
which enforces T{jffl) = T(—SB) in rings with two connecting leads [52, 53, 54] and 
thus no continuous phase shifts can be achieved. In order to break the phase rigidity, it is 
necessary to reduce the device symmetry by attaching additional leads to the ring [54]. 
The addition of leads increases scattering effects in the cross-junctions and requires to 
model the device in a two-dimensional fashion. The effective two-dimensional potential 
at the crossings is sketched in the right panel of Fig. 10. 

Several recently developed recursive Green's function methods principally allow one 
to compute the transmission through Aharonov-Bohm rings [55, 56] but yield the trans- 
mission matrices only for a single Fermi energy. Time-dependent methods based on 
wave-packet dynamics have been implemented for ring structures [57, 58], but have 
the disadvantage that merely the transmission of a certain pulse is detected. The time- 
dependent approach presented in Ref. [51] allows us to obtain the transmission for a 
whole range of energies and is thus ideally suited to describe experiments scanning a 
wide-range of parameters. The theoretical frame- work for the calculation of the macro- 
scopic voltages and currents through the device is the multi-terminal Landauer-Biittiker 
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FIGURE 11. Aharonov-Bohm interferometer and interference pattern, comparison of experimental 
measurements (middle panels) and theoretical calculation (right panels) for two-different voltage and cur- 
rent probe arrangements (upper panels: non-local setup, lower panels local setup). The gray-scale denotes 
the resistivity of the device. Black and white correspond to maxima and minima of the Aharonov-Bohm 
transmission amplitudes and the evolution of the extremal value tracks the phase. Sudden phase jumps 
are indicated by small arrows. The theoretical calculation are performed for a smaller device resulting in 
a larger magnetic field period. The general features of the experiment and the sudden phase-jumps are 
reproduced by the theoretical simulations and are related to multiple reflections within the Aharonov- 
Bohm device. Adapted from Kreisbeck, Kramer (theory) and Buchholz, Fischer, Kunze (experiment) 
[51]. 

formalism [59], which leads to the following expression for the current from channel i 

el' 00 

I i=u dE E k ( -„^(£)i 2 (/(£,jUi,r)-/(£,My,T)), (44) 

fl J CO - / - 

JTi,ni,rij 

where ti ni j nj denotes the transmission amplitude for scattering from the transverse mode 
rij in arm j into the mode rij in arm /. The Fermi functions f(E,jJL, T) = + 
1) _1 characterize the macroscopic reservoirs at the contacts. The semi-infinite lead- 
channels are assumed to be free of imperfections which cause scattering and allow for 
constructing well-defined asymptotic channel eigen-state. A similar concept stands be- 
hind Wigner's and Eisenbud's R-matrix approach, where an artificial boundary is in- 
troduced at the interface between the asymptotic regions and the interaction/scattering 
region. Within the system the probability current is conserved and all equations can be 
derived by considering the probability flux through a closed surface. The R-matrix ap- 
proach describes stationary [60] as well as time-dependent processes [61]. The channel- 
eigenstates extend in principle along the semi-infinite lead. For numerical applications, 
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FIGURE 12. Schematic sketch of the device. The potential is tracked at the contact observations points, 
where particles are injected and removed in accordance with the potential boundary-conditions. 



we construct a wave packet of finite extent by forming a superposition of plane waves 
along the waveguide with a specific transversal mode of the waveguide [46]. Two differ- 
ent measurement setups can be used to probe the oscillatory behavior of the transmis- 
sion. In the local-setup (lower panels in Fig. 11) the time-reversal symmetry is only par- 
tially broken by the finite voltages and the unavoidable device-imperfections, whereas in 
the non-local setup it is possible to adjust the phase by changing the wave-number. The 
middle panels depict experimental results by Buchholz et al. and show that in all cases 
sudden phase-jumps in the signal occur. Our theoretical calculations (right panels) show 
a similar behavior. The phase-jumps can be traced back to multiple reflections along the 
crossings and a resonant coupling of the eigenmodes of the wave-guide and resonant 
states in the crossings [51]. 



INTERACTING MANY-BODY SYSTEMS: DYNAMICAL POINT 

OF VIEW 



Two aspects complicate calculations for interacting many-body systems: the requirement 
of (anti-)symmetrization of the total wavefunction including spin, and the inclusion of 
the Coulomb forces. Exact results for interacting quantum-systems are difficult to ob- 
tain. The case of two-interacting electrons provides already interesting insights into an- 
tisymmetrization effects and is discussed in detail in Ref. [62] using different analytical 
methods, and also numerical and variational approaches. For a larger particle numbers 
(>100), a complete quantum-mechanical description exceeds available computer power. 
Transport calculations face additional difficulties, since in principle the complete elec- 
trical circuit (including the power source) has to be considered. In practice, at a certain 
point we have to make a cut between the microscopic and the macroscopic description. 
The formulation of the correct boundary conditions at this cut presents a formidable 
challenge and an unsolved problem. Most approaches introduce an asymptotic region, 
where the Coulomb-forces are gradually switched off, but it is not clear if this procedure 
captures the nature of a real physical device. In this section we discuss an intriguing 
problem of condensed matter theory, the classical Hall effect [63, 64]. The formulation 
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FIGURE 13. Program flow to obtain the self-consisten Hall potential from a numerical N-body simu- 
lation using a hybrid GPU-CPU scheme. The injection/removal events are timed and spatially distributed 
to maintain an equipotential surface at the contacts. Adapted from Ref. [64]. 



of realistic boundary conditions is already required in a classical many-body theory. The 
computational power for solving classical many -body problems has jumped to new lev- 
els over the last years due to the availability of general purpose Graphics Processing 
Units (GPU). Driven by the ever increasing demand for realistic rendering of computer- 
graphics in games, the processor development of the GPU has overtaken speedwise the 
standard CPUs. Conceptionally, GPUs consists of several hundred stream processing 
units (on an NVIDIA C2050 board 448 units), which work in parallel. The total number 
of floating point operations per second exceeds 1 .2 x 10 12 . GPU programs can be written 
in the Open Compute Language (OpenCL), which allows to use a C++ syntax for the 
main-program, while the GPU routines are coded in a special language. For the Coulom- 
bic problem at hand the fast build-in reciprocal square-root operation of the GPU yields 
another significant speed gain. 
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FIGURE 14. Snapshot of the electron distribution near the source contact. Arrows mark recently 
injected electrons, which are required to maintain the equipotential within the contact region. 



Modelling a Hall device 



Hall nanodevice can be fabricated from AlGaAs/GaAs hetero structures, discussed in 
the previous section. The Hall-bar is schematically sketched in Fig. 12 and consists 
of several interconnected regions. The positively charged donor-layer is shown at the 
bottom and consists of 8094 charges. At a vertical distance of 10 nm the electron 
layer is situated, where electrons can move in a two-dimensional rectangular area of 
length 2.5 /im and width 1.0 fim. At the left-end (source) and right-end (drain) of the 
electron layer, potential observation points are placed. The observation points track the 
time-evolution of the potential due to all the fixed and moveable charges in the device. 
Depending on the observed potentials, electrons are injected or removed. The program 
flow is shown in Fig. 13 and consists of the following steps: Initially we populate the 
device with randomly distributed electrons. Next we calculate the forces acting on each 
electron and update the electron positions by one time-step. The force on the fcth 
electron consists of the Coulombic forces Fp due to all other N e electrons, N c i donors, 
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FIGURE 15. Extremal values of the potential in the source contact, which converges towards the target 
value of Vsource = —5 mV. 



and the velocity dependent Lorentz force F k L , 

F c = f f n-Tk q 2 f r/-r fc 

k 4k£ £ I r z - r k | 3 47t£ £ £J | r, - r k | 3 

= qi k x&, (46) 

where ^ denotes the magnetic field, pointing perpendicular to the device layer. Each 
electron state-vector contains the position r^, and the velocity r k . The equation of 
motions are integrated with Boris' algorithm [65], 

r*(f + A0=- V-^rO + ttt 1 n )r(t) + * (47a) 

rfc(* + Ar)=r)fc(f)+A*rfc(?+A*), (47b) 

where £0/ = j^- denotes the Larmor frequency. Typical time-steps are 1 / 5000 of the 

cyclotron period, which amounts to At = 5 x 10~ 17 s for a magnetic field of SB = 4 T 
and an effective electron mass m = 0.067m e . The observation points at the source and 
drain contacts r SOU rce,/, rdrain,j recor d the local potential averaged over a time-period At RW 

V(Ti) = - l —— f dt' T -— ^ r , (48) 

where the sum runs over all electrons and donors present in the system. The contacts 
serve as interfaces to electron reservoirs and operate under the condition that an equipo- 
tential is required within the contact-pad. The equipotential condition is physically mo- 
tivated by the model of metallic contacts. In practice, we establish an (on average) con- 
stant potential by injection an electron at the observation point, which shows the biggest 
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FIGURE 16. Upper panel: Time-averaged Hall potential between the source (left end) and the drain 
(right end) contact. Lower panel: current density distribution showing the transport pathway and the 
emergence of two "hot-spots" at opposite corners of the device. In the left hot-spot electrons are injected 
into the device. 



positive deviation from the requested contact potential. Similarly, an electron is removed 
at the observation point which shows the biggest negative deviation from the contact 
potential. If none of the conditions is met, the electron propagation continues without 
injection or removal events. The extremal deviations from the prescribed source-contact 
(here V S ource = —5 mV) value are shown in Fig. 15. Convergence is reached after 0.04 ns 
and the complete system approaches a steady-state, where the total injection rate and the 
total removal rate approaches constant values. During the simulation, the state- vector of 
the complete system is stored periodically each 10000 integration steps. After the sim- 
ulation run, a statistical analysis is performed on a regular grid covering the complete 
device area. At each grid point the time-averaged potential V av is obtained, typical av- 
eraging times are in the range of r av = 1 ns, corresponding to 1000 stored state- vectors 
of the system. The time-averaged potential is shown in Fig. 16. The long-time aver- 
age reveals an S-shaped potential, which emerges from in a self-consistent way from the 
many-body calculation. The shape of the potential is a direct consequence of the metallic 
boundary conditions, the particle-interactions, and the specific device geometry. Exper- 
imentally, the same potential has also been observed in the quantum Hall effect (QHE) 
[66, 67]. If we take the observed and theoretically calculated potentials as the mean-field 
potential of the device, we can quantize the system and obtain a theory which contains 
the integer quantum Hall effect [68]. The electric field is strongest in the vicinity of the 
current source contact and leads to a broadening of the local density of states [26] . In the 
time-dependent picture such a broadening signifies a decaying autocorrelation function 



C(t) —> and thus an extended state which guides the electron through the complete 
device to the drain contact. 

CONCLUSION 

The time-dependent approach to complex systems provides an intuitive and dynami- 
cal picture of transport phenomena. The replace stationary state problems by prop- 
agation from initial states, in this way we avoid computationally expensive matrix- 
diagonalization methods, and by the wave-packet method provide results for a wide 
range of energies with a single propagation. For interacting many-body systems, graph- 
ics processing units are powerful tools to devise microscopic models and to test the 
validity of macroscopic equations. 
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